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An Effective Action for Finite Temperature QCD with Fermions 
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Using lattice perturbation theory at finite temperature, we compute for staggered fermions the one-loop 
fermionic corrections to the spatial and temporal plaquette couplings as well as the leading Zn symmetry breaking 
coupling. Numerical and analytical considerations indicate that the finite temperature corrections to the zero- 
temperature calculation of A. Hasenfratz and T. DeGrand are small for small values of k = , but become 
significant for intermediate values of k. The efi'ect of these finite temperature corrections is to ruin the agree- 
ment of the Hasenfratz-DeGrand calculation with Monte Carlo data. We conjecture that the finite temperature 
corrections are suppressed nonperturbatively at low temperatures, resolving this apparent disagreement. The Zn 
symmetry breaking coupling is small; we argue that it may change the order of the transition while having little 
effect on the critical value of /3. 
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1. INTRODUCTION 

It has been known for some time that the effects 
of heavy dynamical fermions can be included in 
Monte Carlo simulations using a hopping param- 
eter expansion of the fermion determinant. This 
is done routinely in perturbative QED where elec- 
tron loops are included in an effective action. Re- 
cently Hasenfratz and DeGrand ||l|[H have per- 
formed a zero-temperature calculation of the shift 
in the lattice gauge coupling constant induced by 
staggered dynamical fermions and applied the re- 
sult to the finite temperature phase transition in 
QCD. Their result for the shift in the critical cou- 
pling, in the form /3P"''° — (3^ + A/J^?, was found 
to hold rather well down to small values of the 
fermion mass, even for Nt ~ 4; this is particu- 
larly surprising since AfSp is calculated using lat- 
tice perturbation theory at zero-temperature. In 
order to understand the effects of finite temper- 
ature, we have calculated the one-loop fermionic 
corrections to the spatial and temporal plaquette 
couplings, as well as the leading Zjv symmetry 
breaking coupling. 
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2. RENORMALIZATION OF /? AT FI- 
NITE TEMPERATURE 

The one-loop finite temperature fermionic cor- 
rection to the 0{A^) term in the gauge field lat- 
tice action is given by 

po 



X A, 



Ap)^[d^}{p) + d\^}{p)-d^^^] (1) 



with 



' sm| ^ I sm 



2 Nt (27r)3 



X S-^{k^pl2)R{K)S-^{k+pl2) (3) 



fen 



(2) 



(4) 



X S-^{k) 

where the vertex functions are given by 

i?(fcp) = cos(fc^) Q(fc^) = -i7p sin(A:^) (5) 



2 



with no sum over /i, and the inverse fermion prop- 
agator is given by 

S{k) ^ ^+ijf,sm{kf,). (6) 

At finite temperature there are two inde- 
pendent tensors of order which are four- 
dimensionally transverse]^. This manifests itself 
as separate renormahzations of the spatial and 
temporal gauge couplings. We find that 
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and 

B{k) = Jj + 4^sin2(fc„). 
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Figure 1. A/? per fermion (spatial and temporal) 
vs. K for various Nt- 



A different result is obtained for A/3f by replac- 
ing 1, 2) in (7) with <i>(fc; 0, 1). Figure 1 com- 
pares A/3 per fermion (spatial and temporal) vs. 
K for TVt = 4 and Nt = 6 with the zero tempera- 
ture result of Hasenfratz and DcGrand. The finite 
temperature values approach the zero tempera- 
ture result from below as a consequence of the 
antiperiodic boundary conditions. Note that they 
ruin the agreement between the zero-temperature 
calculation and the Monte Carlo results for inter- 
mediate to large values of k. 

The connection between the the zero and finite 
temperature result can be understood more phys- 
ically by transforming the sum over Matsubara 
frequencies into a sum over images using the Pois- 
son summation formula for antiperiodic boundary 
conditions. The shift in the spatial plaquette cou- 
pling is then given by 
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with a similar result for A/?t. This form has a 
simple physical interpretation: the integer n la- 
bels the net number of times the fermion wraps 



around the lattice in the temporal direction. Nu- 
merically, the dominant corrections to the zero- 
temperature result come from the first few values 
of n, with the n — \ and n — 2 terms account- 
ing for more than 90% of the finite temperature 
correction for k < 2.0. 
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Figure 2. Diagrams contibuting to finite temper- 
ature corrections for the renormalization of (3. 



The finite temperature corrections result from 
the O(yl^) expansion of image diagrams such as 
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those depicted in Figure 2. The vertical segments 
of the image diagrams are powers of Polyakov 
loops. Given that symmetry is maintained 
in the confined phase of QCD, the contibutions 
of these terms is negligible below /3c- Thus, the 
finite temperature corrections to /3 may be sup- 
pressed nonperturbatively in the confined regime. 
In particular, near (i^ the zero-temperature result 
will hold for A/3. 

3. symmetry breaking in the 
effective action 

There is another set of terms induced by the 
fermion determinant at finite temperature [Q. 
Using the hopping parameter expansion and the 
Poisson summation formula, we find an additional 
contribution to the effective action: 

oo 
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where P{x) denotes a Polyakov loop and the cou- 
plings hn are given by 
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The maximum values of the hn are obtained 
when mp = 0. For A^t = 4 we find h'^'''' = 0.f07, 
^max _ 0.00445, and increasingly smaller values 
for higher order couplings. Including just the hi 
term in the action of an otherwise pure QCD 
Monte Carlo simulation, we observe that this ad- 
ditional source of A/3 is small compared to the 
renormalization of the plaquette couplings dis- 
cussed in the preceeding section. For example, 
even hi = 0.1 at Nt = 4 leads to a shift in /3 per 
fermion of 0.00325. This value of hi corresponds 
to mp = 0.17 which yields a zero-temperature 
predicted shift in f3 per fermion of 0.104. 

The effect of the hi term in the effective action, 
however, is not trivial. We find a first-order phase 
transition for values of hi smaller than approxi- 
mately 0.08. Including zero-temperature plaque- 
tte coupling renormalization, the endpoint of this 



first-order phase transition line maps to the point 
(0.394,4.68) in the (m_F, /3)-plane for the case 
of sixteen degenerate staggered fermion species. 
This is near the endpoint of the first-order phase 
transition observed in simulations with dynami- 
cal staggered fermions ||] . 

A variety of workstations were used to ob- 
tain the results presented in this section. Typ- 
ically 42000 sweeps were made on 10^ x 4 lat- 
tices. Definitive numerics will require substan- 
tially greater computational resources. 

4. CONCLUSIONS 

Finite temperature corrections to the gauge 
coupling renormalization lift the degeneracy of 
the spatial and temporal couplings, but the re- 
sults are in conflict with Monte Carlo data. An 
additional finite temperature shift in /3, due to an 
induced coupling to Polyakov loops, is negligible. 
This coupling does, however, influence the order 
of the transition and may be the most important 
factor in determining the end point of the first- 
order deconfinement phase transition line. As we 
have shown, it is plausible that the success of the 
zero-temperature calculation in determining the 
critical value of /3 is due to a nonperturbative sup- 
pression of finite temperature effects in the low 
temperature regime. Specifically, the small ex- 
pectation value of the Polyakov loops at low tem- 
peratures indicates a suppression of those quark 
paths which account for finite temperature cor- 
rections. 
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